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ABSTRACT 



Context. The region of the Small Magellanic Cloud (SMC) with which this paper is concerned contains the highest concentration 
of IRAS/Spitzer sources, H i emission, and molecular clouds in this neighboring galaxy. However, it has been the target of very few 
studies, despite this evidence of star fomation. 

Aims. We present the first detailed study of the compact H ii region N33 in the SMC by placing it in a wider context of massive star 
formation. Moreover, we show that N33 is a particularly interesting candidate for isolated massive star formation. 
Methods. This analysis is based mainly on optical ESO NTT observations, both imaging and spectroscopy, coupled with other archive 
data, notably Spitzer images (IRAC 3.6, 4.5, 5.8, and S.Oyum) and 2MASS observations. 

Results. We derive a number of physical characteristics of the compact H ii region N33 for the first time. This gas and dust formation 
of 7"4 (2.2 pc) in diameter is powered by a massive star of spectral type 06.5-07 V. The compact H ii region belongs to a rare class 
of Hii regions in the Magellanic Clouds, called high-excitation blobs (HEBs). We show that this Hii region is not related to any star 
cluster. Specifically, we do not find any traces of clustering around N33 on scales larger than 10" (~ 3 pc). On smaller scales, there is 
a marginal stellar concentration, the low density of which, below the 3cr level, does not classify it as a real cluster. We also verify that 
N33 is not a member of any large stellar association. Under these circumstances, N33 is also therefore attractive because it represents 
a remarkable case of isolated massive-star formation in the SMC. Various aspects of the relevance of N33 to the topic of massive-star 
formation in isolation are discussed. 
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1. Introduction 

The Small Magellanic Cloud (SMC) provides a unique 
opportunity for studying star formation in low-metallicity 
environments. Although a neighboring galaxy, the SMC 
is considered to be reminiscent of high-redshift galaxies. 
It has therefore been the subject of numerous surveys and 
studies at various wavelengths. However, as far as massive 
star formation is concerned, detailed study of its individual 
star-forming regions has focused on a relatively small number 
of cases. Apart from the pre-eminent region N66, to which 
a considerable amount of resear ch work has been devoted 
dHev dari-Ma laveri & Seheiil2010i and references therein), only 
sever al star-forming regions as sociated with the Shapley's win g 
JShapley . 1940"), such as N81 d Hevdari-Malav eri et al.1. 1 1999ah 
N88A (iHevdari-Malayeri e t alF 1999b), N83/N84 dLee et al.L 
l2009l) . and N90 dCignoni et al., 2009) have been studied in more 
detail. 

In particular, the star-forming regions populating the lower 
portion of the SMC bar, situated south-west of N66, have 
so far received little attention. The southern end of the SMC 
has the particularity of containing the highest concentration 
of molecular clouds, IRAS/Spitz er sources, and Hi emis- 
sion in this neighboring galaxy (iFukui & Kawamiiral 1201 Ot 



iMizuno et al.','2001';'Stanimirovic et al 



2000 [IT999tlBhtzetal.L 



* Based on observations obtained at the European Southern 
Observatory, La Silla, Chile, Program 69.C-0286(A) and 69.C-0286(B). 



2007; Bolatto et al., 2007; Leroy et alj, I2007t iBot et al.L 12010^ . 
However, higher spatial resolutions and more extended map- 
pings, especially in CO observations, are required to precisely 
associate the molecular clouds of that part with individual 
Hn regions and study the interaction with massive sta rs. Th e 
Hii region LHA115-N33, or in short N33 (HeniH [19561) . 
with J2000.0 equatorial coordinates Q' = 00h49m29s and 5 = - 
73° 26' 34" belongs to this pait of the SMC. Veiy few studies 
have been devoted to it, despite its several interesting character- 
istics, which the present paper aims to highlight. In particular, as 
the present paper will show, this object belongs to a small class 
of Magellanic compact H n regions, called high-ex citation blobs 
(HEBs; for a review see iHevdari-Malaveriet all [2010.) . N33 
appears also under number 138 in the cata log of emission- line 
stars and planetary nebulae compiled by Lindsav ('1961'). In 
contra st, it is not identified in the optical survey of Davies et aH 
(119761). Likewise, N33 bears nu mber 297 in the catalog of 
iMevssonnier & Azzopardil (1 1993b among the list of proven or 
probable SMC compact or small H n r egions. N33 was d etecte d 
as the IRAS source 00477-7343 (iHelou & Walked [19881) . 
and in far-IR as th e source #28 in the ISO 170 fim catalog 
(IWilke et all l2003h . N33 was also included in several radio 
continuum surveys, including the Australia Telesco pe Compact 
Array (ATCA) observa tions at 3 and 6 cm (FiUpovic et al.L 
119981: llndebetouw et al.L[2004l) . 
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This region of the SM C also hosts a planetary n ebula. 
This was first suspected by iHenize & Westerlund' (l96?), who 
nevertheless could n ot con firm its presence. Subsequently, 
iJacobv & De Marcol (l2002h detected a planetary nebula 
candidate at the position (J2000.0) Q; = 00h49m35.0s and 
5 = -73°26' 36"9 (their object number 7). This means that 
the candidate planetary nebula should lie 27" east of the Hii 
region. The reported diam eter of the planetary nebula is 7" 
dJacobv & De MarcH, |2002|) . Although this planetary nebula is 
probably unrelated to the Hn region N33, it has created some 
confusion. In the CDS Simbad database, N33 appears under the 
heading of planetary nebula. It is also notable that among the 1 1 
references given (covering the 1956-2008 period), only two are 
directly concerned with the planetary nebula candidate. 

As a notable point, this object do es not belong to any 
known stellar association of the SMC (Hodg^ Il985h . In the 
context of isolated massive star formation studies, this makes 
N33 particularly interesting. We attempt to show that N33 is 
an isolated Hn region powered by a massive star, which has 
probably formed in isolation. By isolation, we mean not being 
traceable to an origin in any OB association. Observational 
findings^ suggest that massive stars gene rally form in groups 
(e.g., IZi nnecker & Yorke, 2007; Schil bach & Roseii BOOBi: 
iGvaramadze & Bomansii2008j) . There is, however, a statistically 
small percentage of massive stars (~5%) that form in isolation 
dde Wit et all 120051; iParker & Goodwini |2007|) . It is not yet 
understood which physical conditions favor the formation of 
isolated massive stars and how they form. It is pertinent to know 
whether the same continuous process governs the formation 
of clusters of massive stars in large molecular c louds and 
single stars in smaller molecular clouds (see, e.g., lOev et all 
|2004 and below Sect. 4). Should the answer be negative, it 
would be necessary to know the mechanisms responsible for 
the formation of massive stars in isolation. We discuss several 
aspects of this topic. The study of this family of massive stars 
would provide new insights into a clearer understanding of 
massive stars in general, which is not entirely clear. Apart from 
studying global aspects of isolated massive star formation, it is 
necessary to examine each individual case in detail; and N33 
provides a rare opportunity for this research work. 

This paper is organized as follows. Section 2 presents the 
observations, data reduction, and the archive data used (Spitzer 
data, 2MASS data). Section 3 describes our results (overall view, 
physical parameters, extinction, cluster search, stellar content 
and the field population, N33 as a HEB, and Spitzer SED fitting). 
Section 4 presents our discussion, and finally our conclusions are 
summarized in Sect. 5. 



2. Observations and data reduction 

2.1. NTT Imaging 

N33 was observed on 26 and 27 September 2002 using the ESO 
New Technology Telescope (NTT) equipped with the active 
optics and the Superb Seeing Imager (SuSI2; D'Odorico et al., 
[l998;). The detector consisted of two CCD chips, identified as 
ESO #45 and #46. The two resulting frames were automatically 
combined to produce a single FITS file, whereas the space 
between the two chips was "filled" with some overscan columns 
so that the respective geometry of the two chips was approx- 
imately preserved. The gap between the chips corresponds to 
~ 100 true CCD pixels, or ~ 8". The file format was 4288 x 4096 



pixels, and the measured pixel size 0"085 on the sky. Each chip 
of the mosaic covered a field of 5'. 5 x 2'. 7. We refer to the ESO 
manual SuSI2 for more technical information. 

Nebular imaging was carried out using the narrow-band 
filters centered on the emission lines Ha (ESO #884), H/3 
(#881), and [Oiii] (#882) with basic exposures of 300 sec (two 
exposures each on 26 September, six exposures for Ha, and four 
exposures for H/? and [Om] respectively on 27 September). 
The image quality was quite good during the night, being 
represented by a seeing of 0''6. We constructed the line-ratio 
maps Ha/Hy6 and [Oiii]/H/3 from nebular imaging. We also 
took exposures using filters ESO #811 (B), #812 (V), and 
#813 (R) with unit exposure times of 15 sec for B and V and 
10 sec for R, respectively. The exposures for each filter were 
repeated twice using ditherings of 5"-10"for bad pixel rejection. 

PSF-fitting photometry was obtained for all filters using the 
DAOPHOT package under IRAlfl The magnitudes were then 
calibrated using the photometric calibration package photcal. 
To perform this calibration, seven standard stars, belo nging to 
two L andolt photometric groups SA92 and TPhe (iLandoltl 
IT992h were observed at four different airmasses. This led to the 
determination of the photometry coefficients and zero-points. 
Those coefficients are in good agreement with the indicative 
values displayed on the SuSI2 web page. 

The aperture corrections were calculated as follows. Starting 
from one of the frames, we subtracted all stars except the 
ones used for determining the PSF with the daophot.substar 
procedure, using our preliminary DAOPHOT photometry and 
the corresponding PSF. This leads to a frame with only a few 
bright, isolated stars plus residues from the subtraction. We then 
performed both aperture and PSF-fitting photometry on those 
stars, using the same aperture we used for standard stars. The 
comparison led to aperture corrections of 0.02, 0.04, and 0.03 
mag in B, V, and R, respectively. 

During the photometry process, some slight discrepancies 
between the intensity of the frames were found: this effect was 
considered to be the consequence of episodic variations in the 
sky transparency by 7% at most. In order not to introduce a sys- 
tematic underestimation of star magnitudes when averaging the 
frames, we decided to perform photometry on each individual 
frame. 

By cross-correlating the positions of the sources in the vari- 
ous photometry files, we obtained the mean magnitude (average 
of the 2 mag of each filter) and a decent estimator of the uncer- 
tainty in this magnitude (difference between maximum and min- 
imum magnitudes). Finally, the process yielded the photometry 
of 2400 stars in all three filters. The results for the brightest stars 
towards N33 are presented in Table[T] The whole photometry is 
available in electronic form. 

2.2. NTT spectroscopy 

The EMMI spectrograph (iDekker et all [1986) attached to the 
ESO NTT telescope was used on 29 September 2002 to obtain 
several long-slit stellar spectra. The grating was #12 centered on 
4350 A (BLMRD mode) and the detector was a Tektronix CCD 
TK1034 with 1024^ pixels of size 24 jum. The covered wave- 

' http://iraf.noao.fr 
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length range was 38 10-4740 A and the dispersion 38 A mm 
giving FWHM resolutions of 2.70 + 0.10 pixels or 2.48 + 0.13 A 
for a 1"0 slit. At each position, we took three 10 min exposure. 
The instrument response was derived from observations of 
the calibration stars LTT7379, LTT6248, and LTT7987. The 
seeing condition was 0"8 (fwhm). The identifications of the 
stars along the slits (see Fig.[2]for the orientations) were based 
on monitor sketches drawn during the observations. 

Furthermore, EMMI was used on 26 September 2002 to 
obtain nebular spectra with gratings #8 (4550-6650 A) and 
#13 4200-8000) in the REMD mode and with grating #4 
(3650-5350 A) in the BLMD mode. In the REMD mode, the 
detector was CCD #63, MIT/LL, 2048 x 4096 pixels of 15^ /^m^ 
each. Spectra were obtained with the slit set in east-west and 
north-south orientations using a basic exposure time of 300 sec 
repeated several times. The seeing conditions varied between 
0"7. Reduction and extraction of spectra were performed using 
the IRAF software package. Fluxes were derived from the 
extracted spectra with the IRAF task SPLOT. The line fluxes 
were measured by fitting Gaussian profiles to the lines as well 
as by simple pixel integration in some cases. The nebular line 
intensities were co rrected fo r inters tellar reddening using the 
formulae given by iHowarthI (Il983l) for the LMC extinction, 
which is very similar to that of the SMC in the visible. The 
intensities of the main nebular lines are presented in Table [3] 
where F(A) and I{A) represent observed and dereddened line 
intensities. The uncertainties are indicated by the capital letters 
: A < 10%, B=10-20%, C=20-30%, and D > 30%. 



2.3. Archive Spitzer and 2MASS data 

The Spitzer archive data used in this paper come from the 
S^MC project. This is a project to map the star-forming body 
of the SMC with Spitzer in all seven Infrared Array Camera 
(IRAC) and Multiband Imaging Photometer for Spitzer (MIPS) 
bands. We used the IRAC data, obtained in 2005 May, to build 
a composite image of N33 and also obtain photometry. The 
typical PSF of the IRAC i mages in the 3 . 6, 4j , 5.8, and 8.0 
lum bands is l':66 to l':98 jBolatto et all l2007h . The derived 
photometry for N33 in the 3.6, 4.5, 5.8, and 8.0 fim bands are 
11.64, 11.20, 9.61, and 7.91 mag, respectiv ely, using an rntegra- 
tion aperture of 3 pixels, or 3.6" in radius (ICharmandaris et aU 
[200 8). Measurements with either slightly larger or smaller 
apertures do not affect the color results. As for the MIPS fluxes 
(magnitudes) at the 24 and 70 yum bands, they are 0.13 Jy (4.37) 
and 1 .28 Jy (-0.54) respectively. To examine a large field around 
N33 (-400x400 pc), we used the SAGE-SMC observations. 
SAGE-SMC is a Spitzer Legacy program that has mapped 
the entire SMC with IRAC and MIPS. The fuU mosaics are 
available at the Spitzer Science Center homepage, SAGE-SMC 
Data Deliveries. 

We also used the JHK photometry pro- 
vided by the 2MASS point source catalog 
(|http://tdc-www.harvard.edu/catalogs/tmpsc.htm]jl, as pre- 
sented in Table[T] Note that the embedded stars in the H ii region 
(#1, #2, and #3) ai-e not resolved in 2MASS data so the JHK 
photometry of N33-1 corresponds to the whole N33 compact 
H II region. 



3. Results 

3.1. Overall view 

The images taken with the NTT telescope (Sect. 2.1) have a 
whole area of ~5' x5' corresponding to ~ 90 pc x 90 pc for a 
distance of ~60 kpc, or m-M - 18.94 mag (iLanev & StobieL 
Il994l) . They show a starry field marked by very faint diffuse 
nebulosity running in the area. A close-up is presented in Fig. [1] 
which displays the compact Hii region N33 towards the center 
of the image. It has a mean angular radius, {0a. Og)^^^, of 3"7 
corresponding to a radius of 1.1 pc. Broad-band images in B, 
V, and R (Fig.|2]l show the presence of three stars towards the 
H n region (#1, #2, and #3), whose positions and photometry are 
listed in Table[T] The angular separations between these stars are 
2':5 (#1 and #2), 2':5 (#1 and #3), and 3':5 (#2 and #3). We show 
that the central star #1 is the exciting source of the H ii region 
and therefore physically associated with the nebula. As for stars 
#2 and #3, they may be linked to the Hn region or alternatively 
just random field stars. We refer to Section 4 for a discussion. 
Table 1 also contains results for the brightest stars of the field 
close to the Hii region. Star #6 (also known as SMC 013740 
and 2MASS J0 0493037-7326501) is an SMC supergiant K3I 
(iLevesque et al.l, 12006, and below Sect. 3.4). As for star #4, i.e. 
the brightest in the vicinity of N33, it is a Galactic dwarf G5 V 
(see Sect. 3.4). The cross in dicates the positio n of the candidate 
planetary nebula ( Jacobv & De Marcol l2002h . Although these 
authors measure a diameter of 7" for the planetary nebula, i.e. 
similar to N33's, this object does not appear in our image. The 
Spitzer image, a composition of 4.5, 5.8, and 8.0 fim bands 
(Fig.O, shows N33 as a red nebular object with an interesting 
diffuse arc or plume hovering over its eastern and north-eastern 
side. The nature of this feature is not yet clear 

An interesting aspect of the Hii region N33 is its isolated 
character, which we deal with amply in Sects. 3.4 and 4,. As 
a matter of fact, no star clusters and emission nebulae are vis- 
ible in the whole ~ 90x90 pc field o f the NTT images. We 
note however that at the position of N33. lBica & Schmittl ( ll995h 
identify a compact object, which they classify as an emissio n 
nebula, situated in the nebula DEM44 (Davi es et all Il976l) . 
Nevertheless, a positional comparison of N33 with DEM44 (lo- 
cated at 0:48:58.53 -73:25:39.76, J2000) shows that N33 hes 
about 130 pc away from the center of the nebula, and about 
70 pc from its closest edge, considering its dimensions of 7' x 5', 
as provided by Davies et al. ( 1976). As a consequence, there is 
no apparent association of N33 with DEM44. We also exam- 
ined a larger field, ~ 400 x 400 pc, around N33 using the Spitzer 
SAGE-SMC archive data. Similarly, no conspicuous emission 
source is present out to a projected distance of ~ 200 pc north- 
west of N33. At that position, an extended emission object shows 
up, that should correspond to the SMC Hii regions N22/N23. 
Moreover, a minor source is seen south of N22/N23 at a pro- 
jected distance of ~ 160 pc from N33. 

3.2. Physical parameters 

The total Hfi flux of the compact Hii region N33 was de- 
rived using the following procedure. First we calculated the 
relative H/3 flux in an imaginary 1" slit passing through the 
Hy6 image with respect to the total flux emitted by the whole 
Hi: region. This value was then compared with the absolute 
flux obtained from the spectra. The total H/3 flux thus ob- 
tained was F(Hj6) = 1.05 X 10"'^ erg cm"^ s"'. Studies of the 
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extinction in the LMC and the SMC reveal reddening laws 
that are similar to the average Galactic law for the optical 
and near-IR regions (Howarth, 1983; Pre vot et all 119841: 
iBou chet et al.. 1985). Considering the extinction law for the 
LMC ( Howarth, ,1983,) . we computed the reddening corrected 
intensity /(H^S) =4.91 x 10 erg cm"^ s We derived a 
luminosity of 2.2 x 10^^ erg s"', or 550 Lq, for N33 at H(J3). 
This luminosity corresponds to a flux of 5.4 x lO'*^ HyS photons 
s~\ or a Lyman continuum flux of 4.7 x 10**^ photons s ' for 
the star, assuming that the Hii region is ionization bounded. 
The exciting star needed to provide this flux should have an 
effective temperature of ~ 36,000 K or be of spectr al type about 
06.5-07 V, for Galactic metaflicity (Ma rtins et al.L 12005). This 
may however be a lower limit because of photon loss in a 
density-bounded H n region. 

A number of the derived physical parameters of the com- 
pact H II region are summarized in Table|2] The mean angular 
radius of the H ii region, corresponding to the FWHM of cross- 
cuts through the Ha image, is given in Col. 1 . The correspond- 
ing physical radius, obtained using a distance modulus of m - 
M - 18.94 mag (Laney & Stobie, 1994) is presented in Col. 2. 
The dereddened H/3 flux obtained from a reddening coefficient 
of c(Hj6) = 0.67 is given in Col. 3. This reddening coefficient, 
derived from the mean Ha /H/3 ratio of 4.6, corresponds to the 
whole H II region. It is different from the value found from the 
nebular spectrum (Table|3) because, in contrast, the spectrum 
belongs to a particular position and therefore does not cover 
the whole region. The electron temperature calculated from the 
forbidden-Hne ratio [O iii] /l/l 4363/(4959 + 5007), with an un- 
certainty of 4%, is given in Col. 4. The electron density, esti- 
mated from the ratio of the [S ii] doublet /l/l 67 17/6731, is pre- 
sented in Col. 5. It is accurate to ~ 80%. It is well-known that the 
[S ii] lines characterize the low-density peripheral zones of H ii 
regions (see below for corroboration). Column 6 gives the rms 
electron density, >, calculated from the Hf) flux, the radius, 
and the electron temperature, , assuming that the H ii region is 
an ionization-bounded Stromgren sphere. Furthermore, the total 
mass of the ionized gas, calculated from the <ne> with the pre- 
viously noted Stromgren sphere assumption is presented in Col. 
7. The ionization is produced by Lyman continuum photon flux 
given in Col. 8. 

3.3. Extinction 

The average value of the B aimer decrement towards the H ii 
region N33 is about 4.6, corresponding to Ay = 1.5 map , using 
the extinction law for the LMC with R -3.1 (Howarth, 1983). 
The most extincted part of the H ii region is its northeast border, 
where the Ha/Hy6 ratio reaches a value of 6.5 (Ay = 2.5 mag). 
Interestingly, the IR elongated structure (Sect. 3.1), which runs 
from north to south-east (i.e. from the area indicated by the 
upper circle in Fig.|3]l, happens to be adjacent to this higher 
extinction area. The extinction towards star #1 can be derived 
from a second met hod. O-type stars have an intrinsic color of 
B-y=-0.28 mag (Martins & PleH 12006*). This yields a color 
excess of E(B-V) -0.64- mag or a visual extinction of Ay = 
2.0 mag for star #1, in good agreement with the result from the 
Balmer decrement. 

Moreover, the extinction towards N33 was estimated by a 
third method using radio continuum observations. N33 appears 
as the source B0047-7343 in the Parkes radio continuum sur- 
vey at 3 and 6 cm, which had beam-sizes of 2'. 7 and 13'. 8, 



respectively (iFilipovic et al.L Il998h . Higher resolution observa- 
tions of this object in the same wavelength range were obtained 
by fndebetouw et al. (2004), who in their search for ultracom- 
pact and buried H ii regions in the Magellanic Clouds, used the 
Australia Telescope Compact Array (ATCA) in radio continuum 
emission at 3 and 6 cm with synthesized beams of 1"5 and 2" 
respectively. Since the beam- widths are 1 "5 and 2" respectively, 
the radio continuum observations do not sample the entire H ii 
region. To compare the radio continuum fluxes with that of H/?, 
we corrected the HyS flux. We computed these from the ratio of 
the measured Hf) emission flux with different apertures (1"5 and 
2") to the total N33 flux. The resulting extinction. Ay = 1.9 
mag, although comparable with that obtained using the previ- 
ously mentioned methods in the optical range, may be an un- 
derestimate if the flux is not uniformly distributed over the H ii 
region and the smaller radio lobes miss a part of it. 

3.4. Cluster searcli 

We perform a cluster analysis technique on the stellar photomet- 
ric catalog from our SuS12 imaging, to identify any potentially 
important stellar concentration at the surroundings of the 
compact H ii region N33 that would suggest membership of its 
central star to a stellar cluster Our identification method is based 
on star counts in quadrilateral grids for the construct ion of stellar 
density maps in the area of interest (see e.g., Gouliermis et all 
2000). We performed star counts on our complete photometric 
catalog of 2400 sources. This was done by counting the stars in 
grids according to their celestial positions, as they are defined by 
our astrometry. The method is sensitive to three factors, namely 
1) the stellar numbers, 2) the size of each grid element, which 
defines the smallest possibly detectable cluster (the resolution 
element in the density maps), and 3) the size of the available 
field-of-view (FoV). This is due to the isopleths of the stellar 
density maps whose steps are defined to be equal to the standard 
deviation, cr, of the total average stellar density (of stars per grid 
element) of the considered field. We refer to the latter as the 
"background density". We consider any stellar concentration 
revealed by the isopleths, that correspond to stellar density of 
3cr above the background and higher, as statistically important, 
and therefore as a candidate true cluster 

The application of star counts to the whole observed FoV al- 
lows us to identify any possible membership of N33 to any large 
stellar concentration of size on the order of ~ 100 pc. It should 
be noted that the gap between the chips of SuSI2 introduces a 
lack of stars along it in every photometric catalog derived for 
each filter As a consequence, the final combined photometric 
catalog also contains a lack of stars along a 100 pixel thick 
line vertically crossing the middle of SuS12 frame. While this 
introduces a problem into the consistency of the detection, the 
gap is fortunately away from the H ii region, and hence its effect 
on the cluster detection at the immediate vicinity of N33 is 
not important. We constructed the stellar density maps for the 
SuS12 FoV using a grid of 25 x 25 elements, corresponding to 
a physical scale of about 3 pc per grid element, which defines 
the smallest possible stellar concentration that can be revealed 
in the isodensity maps. These maps are constructed for all the 
stars, as well as for selected "blue stars" with B - V < 0.55 
mag, mostly representing the main-sequence stellar population 
and 'red stars' selected to have colors B - V > 0.55 (see 
color-magnitude diagram of Fig.|7]i, which represent the evolved 
stellar population of the general SMC field. In Fig.|4] the stellar 
isodensity maps at the vicinity of N33, i.e. in the close-up field 
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of Figs. [T] and |2] are shown along with the corresponding stellar 
map. 

As seen from the density maps of Fig. |4] all the isopleths 
around N33 appear below the 3cr density threshold. This clearly 
implies that, considering the factors on which the star counts 
method depends, there is no statistically significant stellar 
cluster around N33, of size larger than 3 pc. Our search for 
stellar clustering around N33 is based on the complete stellar 
catalog constructed from our SuSI2 photometry. This catalog, 
according to the evolutionary models, covers stars with masses 
as low as about 2 M© corresponding to late A-type dwarfs. It 
should be noted that the low completeness in faint stars may 
not allow us to have a significant density peak around N33 
simply because we do not see these stars. Nevertheless, a small 
stellar clump of projected size ^ 18"x26"(5 pc x 7 pc) is 
revealed around N33 at the Icr density level and the isodensity 
maps of Fig. ID suggest that it mainly consists of main-sequence 
stars. However, its low density, reaching the 3cr level only 
at its very central peak, does not classify it as a real cluster 
Moreover, from our cluster analysis in the whole observed FoV 
we assess that N33 does not seem to belong to any larger stellar 
concentration on a 100 pc length-scale. A possible reason for 
not identifying any parent stellar structure for N33 is that the 
H II region may indeed be part of a larger stellar aggregate, the 
size of which is much larger than our observed FoV. In this 
case, our search for important density enhancements would 
take place within the system itself, without apprehending its 
existence, thus our average background density would not 
represent the real background, but the high stellar density levels 
of the large conce ntration . Howe ver, p revious searching - by-ey e 
investigations by IHodgd (Il985h and iBica & Schmitf (Il995h . 
on FoVs larger than ours do not detect any large stellar con- 
centration towards this region. The only indi cation of a larg e 
stellar concentration in this region is found by Battinellil (Il99ll) . 
who identified a candidate stellar association, that seems to 
coincide with the larger nebula DEM44. This author applied a 
cluster analysis technique, the path linkage criterion, to previous 
catalogs of OB-type stars to identify candidate OB associations 
in the SMC. However, this method is applied only to the OB 
stars, and therefore is biased towards large-scale "connections" 
of these stars without any supporting information about their 
physical relation or any interaction to each other Moreover, the 
identification of large OB groups based solely on the positions 
of the stars can be challenged, because it has been confirmed 
that young stellar systems may lose a significant fraction of 
their massive stars at the very beginning of their evolution, (e.g., 
iGvaramadze & Bomansl 12008). and therefore OB stars may be 
mistakenly considered as members of a system, while being 
runaways from another Under these circumstances, and because 
DEM44 does not seem to contain N33, the nebula coinciding 
with a candidate OB association cannot be considered proof that 
N33 does belong to this large stellar concentration. 

To assess the clustering behavior of the observed stars in 
the region of N33, we construct the two-point correlation func- 
tion (TPCF), which determines the distance between all possi- 
ble p airs of stars. This method, as applied by e.g. Gom ez et 
(II99 3) to study the distribution of young stars in a Galactic star- 
forming region, considers the excess number of pairs in the ac- 
tual distribution over a random distribution. Here, we apply the 
method as extended by Larson (1995), who introduced the as- 
sessment of stellar clustering in terms of the average surface 
density of companions, 'Ldd), as a function of the projected an- 



gular separation, 6. We measured £c(0) and correlated it to the 
corresponding angular separation for the blue stars in the whole 
observed FoV, as well as for the stars covered by the close- 
up field centered on N33. The constructed TPCFs are shown 
in Fig. |5] drawn with a red line for the close-up field and blue 
for the whole observed field. In both samples, a clear change in 
the correlation 1,^(0) appears at 6 ^ 10", which at the distance 
of the SMC corresponds to about 3 pc. We note that this limit 
is similar to that measured for the small stellar clump seen at 
the position of N33 in the density maps of Fig. |4] For larger 
separations {0 > 10"), the relation 1.^(0) becomes almost flat, 
with oc 0"*'^. This almost-constant density of companions for 
larger separations indicates a random distribution of stars and 
therefore no clustering of stars at scales larger than 3 pc. On the 
other hand, for smaller scales the average density of compan- 
ions correlates linearly with the separation as Ec <^ 0"' ', provid- 
ing clear evidence of self-similar clustering on scales as small 
as the smallest resolved separation, and eventually in the binary 
and multiple systems regime. This implies that stellar clustering 
in the observed FoV seems to occur only at scales smaller than 
about 3 pc, in agreement with the findings of our cluster analy- 
sis above. In conclusion, from our analysis we find that the H ii 
region N33 does not seem to belong to any large stellar concen- 
tration. On smaller scales, we identify a peak in the stellar den- 
sity around N33, with a stellar mass of roughly 100 Mq, which 
nevertheless is very loose and therefore should be considered as 
a random density fluctuation rather than a real cluster 

3.5. Stellar content and the field population 

We obtained two spectra in our program of stellar spectroscopy 
(Sect. 2.2). However in spite of the relatively good seeing 
conditions (0"7 fwhm), extracting uncontaminated spectra is 
not straightforward. The compact Hn region (mean size 7"4) 
has strong emission lines that fill in the absorption lines of 
the embedded stars, in particular those of He i. Nevertheless, 
He II absorption lines at AA42QQ, 4541, and 4686 are certainly 
present, albeit with weak S/N ratio, in the spectra uncorrected 
for nebular lines (Fig.|6). This agrees with the presence of an 
06.5-07 V spectral type, as inferred from the HyS flux measure- 
ment (see Sect. 3.2). The estimated spectroscopic mass of an 
06.5 V star is 29 Mp an d its eff'ective temperature T^tf = 36,800 
K (Mart ins etall l2005i') . We also note that Inde betouw et alj 
(2004) derive a spectral type of 08.5 V for the exciting source 
of N33 from their radio continuum observations at 3 and 6 cm. 
As underlined above, these ATCA observations do not sample 
the entire H ii region. Furthermore, the presence of the He ii 
absorption lines in the spectrum of the exciting star rules out a 
confusion of N33 with a planetary nebula. 

The spectrum of star #4, the brightest object close to the 
Hn region, is pres ented in Fig.[6l Us i ng classification criteria 
for la te-type stars (iJaschek & JaschekL Il987t iGrav & Corballvl 
l2009i') . stai- #4 can be assigned the subtype G5. This spectral 
type is based on the strength of the Can /14227 line and the 
ratios Fe i /14046/H5 and Fe i /14325/Hy. This result agrees with 
the measured colors: B-V^ 0.78, V-R ^0.39,V-K^ 1 .76, 
H - K = 0.07, and J - K = 0.42 mag ( Koornneefl Il983h . In 
terms of luminosity class, it is a dwarf, as suggested by the line 
ratios Srii ^4077/Fe i ^4046 and Yii ^4376/Fei /14383. We also 
note that the Can K and H lines do not have extremely broad 
and damped wings as expected in supergia nts. This star should 
have an absolute magnitude My - +5.1 (iJaschek & JascheR 
Il987h . However, this value is at odds with the observed absolute 
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magnitude of My — -7.31 assuming that it belongs to the SMC. 
This discrepancy indicates that star #4 is a Galactic member. 
Additional evidence of the Galactic membership of this star 
is its relatively low radial velocity. In contrast, star #6, the 
other bright field star, has a radial velocity of 156.4 km s ' 
which confirms its SMC membership. Figure|6] also displays 
the spectrogram of star #6, which is much redder (V - K - 
4.17 mag, see also TablelTJ. We classify this star as an early-K 
type supergiant, i n agreement with its previously reported 
classification, K3 1 ( iLevesque et ali[2006h . 

Figure |7] displays the color-magnitude diagram of the pop- 
ulation of stars in the entire NTT field of 304" x 316" (90 x 93 
pc), for a cut-off" magnitude of V = 20.5. We also overplotted 
four isochrones with ages 18 Myr, 56 Myr, 500 Myr, and 1 
Gyr for metalHcity Z = 0.004 (Lejeune & Schaerer, 2001). A 
fifth isochrone of 14 Myr is calculated with a higher metallicity 
of Z = 0.008. The bulk of the stars are concentrated into two 
groups: an apparent main sequence centered on B - V ~0 mag 
and an evolved population centered on B - V ~0.8 mag. The 
presence of star #6, which is an SMC supergiant of type K3 1 
(explained above), provides a practical constraint in choosing 
the suitable isochrones. Star #6 is appropriately fitted by two 
isochrones, those of 18 Myr and 14 Myr The latter is the 
best fit, but assumes a different metallicity from the other 
isochrones, i.e. a factor of two higher than the average value 
for the SMC. This higher metallicity may represent the effect 
of heavy element enrichment as the K3 star evolves to higher 
luminosities. These 14 Myr and 18 Myr isochrones also closely 
reproduce the properties of the main sequence stars, which in 
their vast majority have an initial mass of roughly < 8 Mq. 

As far as the exciting star of the H ii region N33 is concerned, 
i.e. star #1, it is indicated by a cross on the color-magnitude 
diagram. This star is affected by an extinction of Ay = 2.0 
mag, so it should not belong to the above-mentioned low- and 
intermediate-mass population, which are mostly spread around 
B - V~Q mag. Star #1 appears to be a young massive star 
occurring in a field of old low-mass stars. We also note the 
color-magnitude diagram position of star #2, which apparently 
lies across the Hii region N33. If this star really belongs to the 
56 Myr isochrone, as suggested by the diagram, this rules out 
its association with the H ii region. As for star #3, it seems to be 
consistent with the position of a reddened B-type star, hence its 
association with the H ii region cannot be excluded at this stage. 
The second group of stars on the color-magnitude diagram, 
whose most evolved members gather between B - V colors 0.5 
and 1.5 mag, have a possible age ranging from 500 Myr to 1 
Gyr and a turn-off of ~ 2 Mq. A large number of stars belonging 
to these older isochrones apparently coincide with the lower 
part of the 18 Myr and 14 Myr main sequences. We suggest that 
these various stellar populations may not be locally related to 
star#l. 



an isolated H ii region, the various stellar populations present in 
the color-magnitude diagram should not be physically related to 
it. They are probably distributed along the line of sight over 32 
kpc. We suggest that the Hii region N33 may not be physically 
associated with this main sequence. 

3.6. Chemical abundances 

An example of the nebular spectra of N33 is shown in Fig.|8] 
whereas Table|3] lists the corresponding main lines. The spec- 
trum represents an H ii region, with no He ii lines. The two faint 
lines redward of He i /1 5876 are /1 5888 and A 5890 in emission. 
Both exist in the Orion H n re gion, the fi rst one is unidentified 
and the second is a [Cn] line ^Baldwin et al., 2000). Similarly, 
these two lines are present in the spectrum of the SMC compact 
H II region N66A. The ionic abundances with respect to H^ were 
calculated from nebula r lines using the IRAF task ionic of the 
package NEBULAR (Sha w&DufouR Il995h . The abundance 
values are listed in Table|4] 

To derive the total abundances of a given element, it is 
necessary to estimate the amount of the element in ionization 
states not observed in our spectra. We therefore used a set of 
ionization correction factors (ICFs) to convert into elemental 
abundances. The absence of the He ii line indicates that He^^/H^ 
is negligible. Moreover, we assume that neutral helium is not 
important. Thus we assumed that the total He/H ratio is just 
equal to He+/H+. The total abundance of oxygen was adopted 
to be the sum of and O^^ abundances. The absence of 
Hen recombination lines in our spectra and the similarity 
between the ionization potentials of He^ and O^^ implies 
that the contribution of O^*^ is not significant. To obtain the 
total abundance of nitrogen, we used the usual ICF based 
on t he similarity between the ionization potential of N+ and 
0+ dPeimbert & Costeroi [1969). The N+ abundance does 
not depend strongly on the electron temperature. The largest 
errors come from the uncertainty in the /l/l 6548 and 6584 line 
intensities. Our N result is accurate to within about 30%. As for 
Ne, the only measurable lines in the optical range are those of 
Ne^^ but the amount of Ne^ may be important in the region. We 
adopted the usual expression of the ionization correction factor 
of Ne that a ssumes that the ionization s tructure of Ne is similar 
to that of O (iPeimbert & Costeroi Il969l) . 

The total chemical abundances for N33 are presented in 
Table|5] The most accurately estimated abundances belong to 
He and O, which are accurate to within 15 and 20% respectively. 
Table|5] al so presents the mean a bundance values derived for 
the SMC dRussell & DopitaL[T992 ). The N33 abundances agree 
with the SMC mean values. In other words, N33 is, as expected, 
a low-metallicity ionized nebula. 



The reason is that the SMC is known to have a complex 
structure with considerable line-of-sight depth, as shown by sev- 
eral works; we refer to Westerlund (1990) for a review, as well 
aslMcGee & Newton ( 1981, 1982D. IStavelev-Smith et al.l(ll997h 
andlHatzidimitriou et al.. (2005). In particular. iMathewson et alJ 
(figM l 19881) measured the distances of 161 Cepheids in the 
SMC, using the period-luminosity relation in the infrared. The 
Cepheids were found to extend from 43 to 75 kpc with a max- 
imum concentration at 59 kpc. Since N33 is not part of an OB 
association (see above Sects. 1 and 4) and more importantly it is 



3.7. N33asaHEB 

N33 should belong to the distinct and rare class of Hii 
regions in the Magellanic Clouds called HEBs (see 
iHevdari-Malaveri et al. , 1201 Ol for a review). In contrast to 
the typical Hii regions of the Magellanic Clouds, which are 
extended structures with sizes of several arc minutes corre- 
sponding to physical scales of more than 50 pc and powered 
by a large number of exciting stars, HEBs are relatively dense 



6 



Seller et al.: SMC compact Hn region N33 



Table 1. Positions and photometry of the main stars in the field of SMC N33 ' 



star ID 


a (2000) 


5 (J2000) 


V 


B-V 


V-R 


/ 


H 


K 


spectral type 


N33-1 


00:49:29.25 


-73:26:33.99 


15.92 


0.36 


0.18 


13.94 


13.31 


12.96 


06.5-07V 


N33-2 


00:49:29.62 


-73:26:35.70 


16.39 


1.07 


0.61 










N33-3 


00:49:28.84 


-73:26:35.28 


17.50 


0.28 


0.10 










N33-4 


00:49:24.60 


-73:26:45.51 


11.65 


0.78 


0.39 


10.31 


9.96 


9.89 


G5 V (Galactic) 


N33-5 


00:49:41.75 


-73:26:41.83 


13.04 


0.77 


0.40 


11.72 


11.38 


11.33 




N33-6* 


00:49:30.58 


-73:26:50.87 


13.56 


1.76 


0.95 


10.48 


9.62 


9.39 


K3I 


N33-7 


00:49:36.86 


-73:27:54.49 


13.89 


0.59 


0.36 


12.55 


12.19 


12.09 




N33-8 


00:49:14.22 


-73:26:15.80 


14.71 


1.48 


0.33 


12.10 


11.33 


11.18 




N33-9 


00:49:23.10 


-73:25:53.87 


15.01 


0.36 


0.20 


14.13 


13.86 


13.82 




N33-10 


00:49:43.17 


-73:26:37.02 


15.03 


0.18 


0.11 


14.47 


14.30 


14.28 




N33-11 


00:49:18.73 


-73:27:31.13 


15.23 


1.27 


0.66 


12.96 


12.32 


12.21 




N33-12 


00:49:39.42 


-73:27:17.54 


15.27 


-0.19 


-0.13 


15.75 


15.66 


15.45 




N33-13 


00:49:39.61 


-73:26:13.45 


16.10 


1.37 


0.72 


13.67 


12.93 


12.83 





t The BVR photometry results from NTT observations while the JHK measures come from the 2MASS catalog. 
t Also named SMC 013740. 



Table 2. Some physical parameters of the compact H ii region SMC N33 



d 


r 


Km' 


Te 






^gas 


Nl 


n 


(pc) 


erg s"' cm"^ 


(K) 


cm"^^ 


cm"^ 


(Me) 


ph s"' 






X 10-'2 










xlO'*'* 


3.7 


1.1 


4.91 


12540 


450 


380 


64 


4.7 



t Corrected for reddening with c(Hy8 ) = 0.67. 
:|: Estimated from the [S ii] ratio. 



and small regions of ~5" to 10" in diameter in the optical, 
corresponding to ~ 1.5 to 3.0 pc and excited by a much smaller 
number of massive stars. Their excitation, as derived from their 
[O m] (AA 4959 + 5007) / Hy8 ratio, is generally larger than that of 
ordinary Magellanic Cloud H n regions. For a fixed metallicity, 
the [Om]/Hy6 ratio increases with the effective temperature of 
the exciting star, as well as with the gas density in homogeneous 
photoionization models. 

N33 qualifies for membership of this category. It is compact, 
with a mean angular radius of 3"7, or 1.1 pc. The compactness 
impUes a young age, since H ii regions disperse with time and 
become diff'use. It has a mean [Oiii] /HyS ratio of 5.9, which is 
higher than the mean value for ordinary extende d H ii regions in 
the Magellanic Clouds dHevdari-Malaveri et all l2010i and ref- 
erences therein). The electron temperature amounts to 12,540 
K, which corresponds to excitation by a massive star at least as 
hot as an 06.5-07 type at SMC metallicity. The high-excitation 
auroral transition of doubly ionized oxygen, [Om]/l4363, is 
quite strong, as shown in Figs.|6]& 8. Moreover, the estimated 
mean extinction over the H ii region, as derived in the optical, 
is Av = 1.5 mag. This extinction peaks at 2.5 mag towards 
the northeast border, where the nebula is probably in contact 
with a molecular cloud of unknown characteristics. This as- 
sociation with a molecular cloud is another characteristic of 
HEBs. Moreover, in the plot comparing Magellanic Clouds com- 
pact Hn regions among them, N33 appears to be one of the 
most excited and brightest members of this class in the SMC 
dMevnadier & Hevdari-Malaveri. ,2007,) . 

3.8. Spitzer SED fitting 

We used our Spitzer photometry to construct the mid-IR SED 
of N33 and assess the behavior of the circumstellar dust of the 
nebula. We fitted this SED with the library of YSO models by 



iRobitaille et al.l (l2006h using the onUn e SED fitting tool, pro- 
vided by these authors ("Robitaill e et al.L |2007l^. Figure|9] shows 
the ten best-fit models to this SED. Only three of these mod- 
els include both circumstellar envelopes and disks. The majority 
of the best-fit YSO models contain only one envelope around 
the central stellar source of N33 with an outer radius varying 
between 3x10"* and 10^ AU. The derived total luminosity of 
the nebula according to these models varies between 5x10^ 
to IO'^Lq. This agrees with the luminosity of ~5 x 10'* L©, de- 
rived from IRAS fluxes considering that the IRAS SED peaks 
at ~ 60 fim. Unfortunately, the blue part of the complete SED 
of N33, which includes the stellar source, could not be prop- 
erly fitted. There are four models available in the library of 
IRobitaille et al. I (l2006h that include stellar sources with T^^g, cor- 
responding to the spectral type of the central star of N33 and a 
mass of about 25 Mq. However, their fit to the blue part of our 
SED was not successful because there were few available mea- 
surements in the optical. In all cases, although SED fitting is an 
interesting exercise, its results should not be over-interpreted for 
the simple reason that N33 is not a YSO. 

4. Discussion 

N33 is quite noteworthy in that it is a isolated Hii region in 
this area of the SMC. In contrast to typical Hii regions in the 
Magellanic Clouds, which have neighboring Hii regions, N33 
differs from any typical emission nebula. Thus it constitutes 
an interesting case for studies regarding the physical condi- 
tions in which massive stars form in isolation. Observational 
findings suggest that massive star formation is a collective 
procesSj_i.e. that massive stars form in cluster environments 
(e.g.. IzTnnecker & Yo rke. 2007; Schilbach & Roser, 20081 
iGvaramadze & Bomansl l2008h and that the mass of the most 
massive star in a cluster is correlated with the mass of the 

- and available at http://caravan.astro.wisc.edu/protostars. 
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Table 3. Nebular line intensities of the SMC compact H ii region N33 



A ( A) 


Men. 


F(A) 


1(A) 


Accuracy 


3727,29 


[On] 


11. 1 


108.8 


A 


3797 


HIO 


2.4 


3.3 


B 


3835 


H9 


3.0 


4.1 


B 


3869 


[Ne ra] 


15.2 


20.3 


A 


3889,90 


Hei + H8 


8.6 


11.4 


B 


3968,70 


[Ne ra] + He 


11.8 


15.2 


A 


4101 


m 


14.0 


17.2 


A 


4340 


Hy 


29.1 


33.3 


A 


4363 


[Oni] 


4.8 


5.5 


B 


4471 


He I 


3.4 


3.7 


C 


4861 


HyS 


100.0 


100.0 


A 


4959 


[Oni] 


148.3 


145.1 


A 


5007 


[Oni] 


444.8 


430.8 


A 


5876 


He I 


13.8 


11.5 


B 


6300 


[Oi] 


4.4 


3.4 


C 


6312 


[Sni] 


2.3 


1.8 


C 


6363 


[Oi] 


1.3 


1.0 


D 


6548 


[Nn] 


4.1 


3.1 


C 


6563 


Ho- 


380.0 


286.0 


A 


6584 


[Nn] 


10.3 


7.8 


B 


6678 


He I 


4.1 


3.1 


B 


6716 


[Sn] 


10.1 


7.5 


B 


6731 


[Sn] 


9.3 


6.9 


B 


7065 


He I 


4.3 


3.1 


C 


7135 


[ Ar ni] 


11.9 


8.4 


B 


7236 


[ Ar iv] 


1.0 


0.7 


D 


7323 


[On] 


10.1 


7.0 


C 


7751 


[ Ar ni] 


3.7 


2.4 


C 


c(I^) = 0.40 



Table 4. Ionic abundances of SMC N33 



Ion N33 

HeVH+ 0.087 

0+/H+ (X 10^) 3.65 

0++/H+ (X 10') 7.45 

N+/H+ (X 10^) 0.92 

Ne++/H+ (X 10*) 8.89 

S+/H+(xlO') 2.10 

Ar++/H+ (X 10') 4.73 



Table 5. Elemental abundances in SMC N33 



Element 


N33 


mean SMC^ 


He/H 


0.087 


0.081 


O/H (X 10") 


1.11 


1.07 


N/H (X 10*) 


2.76 


4.27 


Ne/H (X 10') 


1.34 


1.86 



t Russell fePopita (1992) 



cluster itself dElmegreenl l2000t IWeidner & Kroupat l2006h . 
Nevertheless, other observational results imply that massive 
stars do not necessarily form in clusters but can be formed 
either as isolated s tars or in very sma ll groups. According to a 
statistical study bv lde Wit et all (12005 ). nearly 95% of Galactic 
O population is located in clusters or OB associations. This 
means that a small percentage, about 5%, of high mass stars may 
form in isolation, in very good agreement with the finding by 
iParker & GoodwinI (l2007h . Even in the Galactic center region, 
a number of Wolf-Rayet s tars and O supergiants might have 
formed in isolation dMauerhan et al.L 12010.) . ,Lamb et al.l (1201 Ot) 



observed eight apparently isolated O stars in the SMC. Among 
the six non-runaway cases, there is no evidence of clustering in 
three of them, which remain interesting candidates for isolated 
massive-star formation. By isolation, we mean not clearly 
originating in an OB association. This definition therefore 
excludes runaway massive stars, which are thought to result 
from either dynamical interaction in massive dense clusters, 
or via a kick from a sup ernova explosion in a binary system 
(e.g., iGvaramadze et akf |2009; Pflamm-Altenburg & Kroup^ 
1201 Ol and references therein). Parker & GoodwinI (12007) give 
a different definition for an isolated massive star: an O-type 
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Fig. 1. A composite three-color image of the area of the SMC containing the Hii region N33. The Hii region the bright object is 
situated above the image center at coordinates a - 00h49m29s and 6 = -73°26' 34" See Fig.|2]for the identification of other objects. 
The image, taken with the ESO NTT/SuSI2, results from the coaddition of narrow-band filters Ha (red), [O iii] (green), and H/3 
(blue). Field size 156"x 147" (~ 46 x43 pc). This is a close up of an original image covering a field of 304" x 3 16" corresponding 
to 90 X 93 pc. North is up and east to the left. 



star belonging to a cluster whose total mass is <100Mq and 
moreover is devoid of B-type stars (10 < M/Mq< 17.5). In 
the LMC, s everal massive runaway candidates have been 
reported (e.g. Hevdari-Malaveri et all l2000l lEvans et all 1201 Ot 
[Cjvaramadze et al., 2010a). Similarly, in the SMC a dozen OB 
runaway candidates were reported via detection of their bow 
shocks using Spitzer data (iGvaramadze et alll2010bh . 

In any case, the massive star(s) powering N33 cannot be 
runaway(s) for three main reasons: 1) N33 is not linked to 
any OB association; 2) it seems impossible that a massive star 
carries its H ii region during the ejection, although ejection into 
a molecular cloud cannot be a priori excluded; and 3) were N33 
to contain more than one (probably exciting) stars, this situation 
would make the ejection hypothesis still less probable. On the 



other hand, the possibility that a massive star ejected into a 
molecular cloud creates an H ii region around itself has never 
been encountered. 

To obtain a clear picture of how scattered the cold inter- 
stellar material is in the general area of N33, we investigate 
the distribution of molecular clouds and known young stellar 
objects (YSOs) in the vicinity (~ 200 pc) of N33. CO surveys 
of the SMC reveal a concentration of molecular clouds in 
the southern end of this galax v (iFukui & K awamura!. 1201 Ol: 
Mizuno et al., 2001; Stanimirov ic et al.L 2000, 1999; Blit z et all 
2007: iBolatto et aU i2007l: iLerov et al.ri2007c iBot et al.1. l2010h . 
However, the available CO observations do not detect the 
N33 Hn region, maybe because of a lack of sensitivity 
or resolution or for other reasons. The nearest molecular 
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1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 — 

45.0 40.D 35.0 0:49:30,0 Z5.D ZO.O 15.0 

Right Ascension (JZOOO) 

Fig. 2. A broad-band image of the SMC N33 region obtained through filter R (ESO #813). Same field size and orientation as in 
Fig-El The arrows 1 and 11 show the directions of the spectrograph slit. Seeing 0"6 FWHM. The Hn region apparently contains 
three stars (numbered #1, #2, and #3). The brightest stars of the field are also numbered. The cross indicates the position of the 
candidate planetary nebula reported bv lJacobv & De Marcol(l2002h . 



cloud to N33 is cen t ered a t a projected distance of ~230 pc 
(iFukui & KawamuraL |2010", and references therein). As for 
YSOs, we used the catalog of candidates for these objects 
derive d with Spitzer from the S"^MC survey dBolatto et al.L 
12007). In this area, we found five candidate YSOs, three of 
them being quite remote from N33 at distances of between 
about 48 and 80 pc, hence we do not consider them to be related 
to N33. On the other hand, we note that the two remaining 
candidate YSOs, named S3MCJ0 04930.12-732623.4 2 and 
S3MCJ004929.07-732658.56 (see Bola tto et all l2007l Table 
4), are projected much closer to N33 at distances of about 5 
and 7 pc, respectively (Fig.[3]l. The nature of these red sources 
as true YSOs still remains to be verified, and the possibility 
of a connection between N33 and these objects certainly 
requires further investigation, preferably with mid- and near-lR 
spectroscopy of these sources. We also note that the S^MC team 



did not identify N33 as a candidate YSO. 

It has been suggested that isolated field OB stars originate 
in a different m ode of star f ormati on from their cou nterparts 
in associations. iMassev I (120021) and iKroupa & Weidneil (liool) 
make the case that the field star initial mass function (IMF) 
might differ from that in OB associations/clusters. More explic- 
itly. Galactic-field IMFs for early-type stars cannot, under any 
circumstances, be a Salpeter power-law, i.e. a = 2.35 , but must 
have a steeper exponent, a > 2.8 ( Krou pa & WeidnerL l2003h 
Similarly, special physical conditions have been called up to 
account for the formation of very high-mass isolat ed stars in the 
bulge of the interacting Whirlpool galaxy, M51 (Lamer s et all 
2002). According to these authors, the formation of the bulge 
stars in M51 seems to be in line with theoretical predictions 
that isolated massive-star formation might take place in clouds 



10 



Seller et al.: SMC compact Hn region N33 




Fig. 3. A composite image of the SMC N33 region obtained with Spitzer IRAC data. The object at the center is the H ii region N33. 
The two other bright sources are stars #4 and #6 (see Fig.|2). The 4.5 fj.m band is represented in blue, the 5.8 fim band in yellow, and 
the 8.0 ;um band in red. The compact H ii region N33 is the red nebular object at center Note the diffuse emission curl over the main 
body of the H ii region running from north (inside the upper circle) to south-east. The circle s enclose th e young stellar object (YSO) 
candidates S3MCJ004930. 12-732623.42 and S3MCJ004929.07 -732658.56 reported bv lBolatto et al.i (20Q7). However, the true 
nature of these objects and their relation with N33 are not clear; see the text (Sect. 4). Field size 280"x 160" (~ 83 x 47 pc). North 
is up and east to the left. 



where H2, [Oi]63yum, and [Cii]158yum are the dominant 
coolants (Norman & Spaans, 1997; Mihos et al., 1999). These 
conditions are expected to occur in regions of low CO and 
dust contents because of the low metallicity, where the optical 
depth is Ay < 1 mag and t he presence of a h ot source can 
dissociate the CO molecules (iLamers et al l l2002h . The isolated 
massive stars might therefore be the result of these physical 
conditions and the luminous nucleus (situated at a distance 
between 90 and 270 pc). For comparison, the metallicity of 
SMC N33 is about a factor of 10 lower than solar. However, no 
major CO-dissociating hot source is pres ent. The p re-eminent 
massive star cluster N66/NGC 346 (Hevd ari-Malaveri _& Selier„ 
I2OIOI and references therein) lies some 1400 pc away. Could 
this source have a similar effect as the M51 nucleus? Another 
similarity with M5 1 is that the SMC is also an interacting galaxy. 

We note that in contrast to the above-mentioned proposal 
for t he two diff e rent f ormation modes for cluster and isolated 
stars, lOev et all (l2004l) find that the formation of field massive 
stars is a continuous process in associations. More specifically, 
the field stars do not originate from any different star-formation 
mode. This conclusion is based on an empirical census of 
uniformly-selected massive-star candidates distributed all-over 
the SMC in more than 100 star clusters containing at least 3 
to about 80 stars. Moreover, according to their Monte Carlo 
simulations, single, field OB stars are usually the most massive 
member of a group of smaller stars. This result is consistent 
with the idea that there is both a universal IMF and universal 



N^^ clustering law, which extends to A^, = 1 . Jointly, these laws 
imply that the fraction of field OB stars typically ranges from 
about 35% to 7% for most astrophysical situations. 

iKauffmann & Pillail (l20Toh suggest that massive stars might 
form only in molecular clouds possessing a minimum mass 
depending on the cloud size: M(r) > 870 Mq (r/pc)'-^^ Thus, 
massive-star formation requires a large mass being concentrated 
in a relatively small volume. Clouds below this threshold might 
still form stars and clusters of up to intermediate mass (such as 
the Galactic Taurus, Ophiuchus, and Perseus molecular clouds). 
One can infer from this finding that molecular clouds with 
a mass slightly above this limit can form one or a couple of 
massive stars. N33 may represent such a case, but the problem 
is not at all solved, since it is then unclear why in this area of the 
SMC only one cloud happens to be massive enough to give rise 
to massive stars. In other words, we have yet to ascertain which 
physical conditions inhibit the aggregation of smaller clumps 
into massive clouds. 

In any case, the isolated massive star(s) powering N33 
cannot be a result of t he competitive accretion process 
(iBonnellet all 119971 12001h . According to this scenario, devel- 
oping protostars in their natal molecular clouds compete with 
each other to gather mass. The protostars accrete mass with a 
rate that depends on their location within the protocluster. They 
use the same reservoir of gas to grow. Hence, the protostars 
nearest the center, where the potential well is deep and gas 
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Fig. 4. Isodensity contour maps constructed to search for stellar clustering around the compact SMC H ii region N33. In the top-left 
panel, the stellar map for the close-up field of view centered on N33 (position 0,0). The stellar density maps refer to the same 
close-up field of N33. They are shown as they are constructed for all observed stars (green map, top-right panel), for blue (main 
sequence) stars (blue map, bottom-left panel), and red evolved stars (red contour map, bottom-right panel). Isopleths are plotted 
starting at the average background stellar density and in steps of cr above this level. No stellar concentration with density > 3cr 
above the background is identified in the maps. Only the central part of the small stellar clump encircling N33 in the green and blue 
maps reaches the 3cr density level. 



densities are higher, have the highest accretion rates. This model 
is synonymous with collective massive-star formation. 



There is no doubt that star #1, the exciting source of the H n 
region, is a massive star (of derived spectral type 06.5-07 V). 
We also note that this star constitutes the main object of the 
isolated massive-star formation debate presented in this paper 
In contrast, the nature of the two neighboring stars #2 {V - 
16.39, B-V^ 1.07 mag) and #3 (V = 17.50, B - V = 0.28 mag) 
is unclear. The presence of these stars raises several questions: 



Are they massive? Are they physically related to the compact 
H II region and in particular to the massive star #1 ? The present 
study cannot provide firm responses to these questions. We have 
found that if there is a clustering in the region this occurs at 
a scale of about 3 pc and smaller, which actually covers the 
1 pc separation between stars #1, #2, and #3. We note, however, 
that the color of star #2 (B - V = 1.07 mag) does not match 
that expected for a massive star assuming that stars #1 and #2 
are more or less equally affected by reddening. As for star #3, 
its absolute magnitude and color are compatible with a B-type 
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Fig. 5. Average surface density of stellar companions as a function of projected angular separation from each blue star in the whole 
observed area (blue solid line) and in the region of the vicinity of the Hii region N33 (red dashed line). Error bars are indicated. 
From both curves, one can distinguish a break in the correlation at almost the same separation of - 10", i.e., 3 pc. This scale 
separates any large stellar clustering from the small binary/multiple stellar clumps. 



star of mass about 15 Mq. Further high-resolution observations, 
both imaging and spectroscopy, are necessary to probe the exact 
nature of this star and its relation to the H ii region. 

HEBs are usually located adjacent to ordinary giant Hn 
regions or seen lying across them. This implies that their 
formation is the consequence of triggering by a previous 
gener a tion of massive st a rs in th e complex (Elmegreen & Ladal 
119771: IWhitworth etsl[ Il994t iDeharveng et al.. 20091 and 
references therein). However, this is not the case for N33 
whose massive star(s) has(have) formed in isolation. Another 
example of an SMC HEB that has formed in isolation is 
N81, which was studied using HST imaging and spectrosc opy 
dHevdari-Malaveri et al.L Il999ai [2002: Martins et al.L l2004t) . It 
seems, however, that N33 is a comparatively more attractive 
candidate for isolated massive-star formation. N8 1 is an isolated 
Hii region in Shapley's wing. Nevertheless, it is certainly pow- 
ered by at least two main O type stars, which are apparently part 
of a small c luster of massive stars d e tected w ith the HST spatial 
resolution (iHevdari-Malaveri et al.L Il999al) . The massive star 
at the origin of N33 has most certainly formed in an isolated, 
small molecular cloud, the mass and physical characteristics of 
which are not known. 



It will be interesting to derive the IMF for such small 
massive star clusters formed in isolation, and compare it with 
that of massive stars formed in OB associations. We emphasize 
that individual studies of a variety of cases is essential in 
parallel with global, statistical investigations. These small 
clusters provide relatively simpler situations for studying the 
IMF because they involve a smaller number of physical factors 
because of their isolation. In the case of N33, it seems that the 
initial mass of the exciting star is ~40Mq. It is unsurprising 
that the highest masses in small clusters be lower than those in 
OB associations. Indeed small clusters originate in molecular 
clouds with smaller masses. To advance the study of the mass 
distribution in small clusters, high resolution techniques in the 
optical and IR, both imaging and spectroscopy, environmental 
as well as stellar, are required. 



5. Concluding remarks 

This paper has presented the first detailed study of SMC N33 
using imaging and spectroscopy in the optical obtained at the 
ESO NTT, as well as from the Spitzer and 2MASS data archives. 
We have derived a number of physical characteristics of N33 
and ascertained its powering source. This compact H ii region of 
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Fig. 6. Spectrograms of three stars observed towards SMC N33. Star # 1 is the exciting source of the compact H ii region N33. This 
is a raw spectrum, i.e. it is not corrected for nebular emission lines. Note the He ii absorption lines indicating a hot massive star (see 
Sect. 3.5). Stars #4 and #6, lying in the field of N33, are classified as G5 V (Galactic) and K3 1 (SMC) respectively. 



~ 7"4 (2.2 pc) in diameter belongs to a small class of HEBs in 
the MCs. In contrast to other members of this class, N33 is not 
associated with any OB association. This object is excited by 
a massive star of type 06.5-07 V, and represents an interesting 
case of isolated massive-star formation in the SMC. 
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